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Abstract 

We prove the global in time existence of spherically symmetric solutions to an 
initial-boundary value problem for a system of partial differential equations, which 
consists of the equations of linear elasticity and a nonlinear, non-uniformly parabolic 
equation of second order. The problem models the behavior in time of materials 
in which martensitic phase transitions, driven by configurational forces, take place, 
and can be considered to be a regularization of the corresponding sharp interface 
model. By assuming that the solutions are spherically symmetric, we reduce the 
original multidimensional problem to the one in one space dimension, then prove the 
existence of spherically symmetric solutions. Our proof is valid due to the essential 
feature that the reduced problem is one space dimensional. 

1 Introduction 

Many inhomogeneous systems can be characterized by domains of different phases sep- 
arated by a distinct interface. When driven out of equilibrium, their dynamics result in 
the evolution of those interfaces, and the systems might develop into structures (com- 
positional and structural inhomogeneities) with characteristic length scales at the nano-, 
micro- or meso-scale. To a large extent, the material properties of such systems are de- 
termined by those structures of small-scale. Thus it is important to understand precisely 
the mechanisms that drive the evolution of those structures. Materials microstructures 
may consist of spatially distributed phases of different compositions and/or crystal struc- 
tures, grains of different orientations, domains of different structural variants, domains 
of different electrical or magnetic polarizations, and structural defects. These structural 
features usually have an intermediate mesoscopic length scale in the range of nanometers 
to microns. The size, shape, and spatial arrangement of the local structural features in 
a microstructure play a critical role in determining the physical properties of a mate- 
rial. Because of the complex and nonlinear nature of microstructure evolution, numerical 
approaches are often employed. For more details, see e.g. [T3l [F4"l [22]. 
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In this article we are interested in a model for the evolution, driven by configurational 
forces, of microstructures in elastically deformable solids. There are two main types of 
modeling for the evolution of microstructures. In the conventional approach, the regions 
separating the domains are treated as mathematically sharp interfaces. The local interfa- 
cial velocity is then determined as part of the boundary conditions, or is calculated from 
the driving force for interface motion and the interfacial mobility. This approach requires 
the explicit tracking of the interface positions. Such an interface-tracking approach can 
be successful in one-dimensional systems, however it will be impractical for complicated 
three-dimensional microstructures. Therefore, during the past decades, another approach 
has been invented, namely, the phase-field approach in which the interface is not of zero 
thickness, instead an interfacial region with thickness of certain order of a small regular- 
ization parameter. Though it is still a young discipline in condensed matter physics, this 
approach has emerged to be one of the most powerful methods for modeling the evolu- 
tion of microstructures. It can be traced back the theory of diffuse-interface description, 
which is developed, independently, more than a century ago by van der Waals [26] and 
some half century ago by Cahn and Hilliard |llj . 

The two well-known models for temporal evolution of microstructures are the Cahn- 
Hilliard/ Allen-Calm equations corresponding, respectively, to the case that the order 
parameter is conserved and not conserved. These phase field models describe microstruc- 
ture phenomena at the mesoscale (see e.g. [22] ) . and one suitable limit of it may be the 
corresponding sharp- or thin-interface descriptions. In this article we study a model for 
the behavior in time of materials with diffusionless phase transitions. The model has 
diffusive interfaces and consists of the partial differential equations of linear elasticity 
coupled to a quasilinear, non-uniformly parabolic equation of second order that differs 
from the Allen-Cahn equation (the Cahn-Hilliard equation in the case that the order 
parameter is conserved) by a gradient term. It is derived in [21 3] from a sharp interface 
model for diffusionless phase transitions and can be considered to be a regularization of 
that model. To verify the validity of the new model, mathematical analysis has been car- 
ried out for the existence/regularity of weak solutions to initial boundary value problems 
in one space dimension, [31 [5j El EZ1 [28], the motion of interfaces [6], and the existence 
of traveling waves [19]. In the present article, the existence of spherically symmetric 
solutions to an initial-boundary value problem will be studied. We first formulate this 
initial-boundary value problem in the three-dimensional case, then reduce it, by assuming 
that the solution is spherically symmetric, to the one-dimensional case. The existence of 
weak solutions to this one dimensional problem is proved. 

Let !] C K 3 be an open set. It represents the material points of a solid body. The 
different phases are characterized by the order parameter S(t,x) € K. A value of S(t,x) 
close to zero indicates that the material is in the matrix phase at the point x € 17 at 
time t, a value close to one indicates that the material is in the second phase. The other 
unknowns are the displacement u(t,x) £ M 3 of the material point x at time t and the 
Cauchy stress tensor T(t, x) 6 S 3 , where <S 3 denotes the set of symmetric 3 x 3-matrices. 
The unknowns must satisfy the quasi-static equations 



div x T(t, x) 
T(t,x) 




b(t,x) 



(1.1) 
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for (t, x) G (0, oo) x Q. The boundary and initial conditions are 

u(t,x) = j(t,x), S{t,x) = 0, (t, x) £ [0, oo) x dQ, 1 (1.4) 

S(0,x) = S (x), ie!l. (1.5) 

Here V x u denotes the 3x3-matrix of first order derivatives of u, the deformation gradient, 
(S7 x u) T denotes the transposed matrix and 

e{V x u) = - [V x u + (V x uf) 

is the strain tensor, e € <S 3 is a given matrix, the misfit strain, and D : S s S 3 is the 
elasticity tensor, a linear, symmetric, positive definite mapping. In the free energy 

i(>(e, S) = \ {D(e - sS)) ■ (e - eS) + $(S), (1.6) 

we assume that ip G C 2 {R, [0,oo)), choose 4> as a double well potential with minima at 
S = and S = 1. is t ne partial derivative. The scalar product of two matrices ^4 and 
B is denoted by A ■ B = ^ dijbij. c > is a constant and is a small positive constant. 
Given are the volume force b : [0, oo) x fi -> M 3 and the data 7 : [0, 00) x <9f2 — > R 3 , 
S : n K. 

This completes the formulation of the initial-boundary value problem. Equations 
(jl.ip and (ll.2p differ from the system of linear elasticity only by the term eS. The 
evolution equation (|1.3p for the order parameter S is non-uniformly parabolic because of 
the term uA.S\V x S\. Since this initial-boundary value problem is derived from a sharp 
interface model, to verify that it is indeed a diffusive interface model regularizing the 
sharp interface model, it must be shown that the equations (jl.ip - (jl.5p with positive v 
have solutions which exist globally in time, and that these solutions tend to solutions of 
the sharp interface model for v — > 0. This would also be a method to prove existence of 
solutions to the original sharp interface model. 

We only contribute to the first part of this program in this work and show that there 
exist some special solutions to the initial-boundary value problem that is essentially in one 
space dimension. Up to now we still can't solve the following problem: either solutions 
in three space dimensions exist or these solutions converge to a solution of the sharp 
interface model for v — > 0. We shall see later that the existence result of spherically 
symmetric solutions is of interest because the problem has a stronger nonlinear term 
(compared with the problem by assuming all unknowns depend on only one component 
of space variable x which is studied in [3]), despite it is essentially one space dimensional. 

Related to our investigations is the model for diffusion dominated phase transforma- 
tions obtained by coupling the elasticity equations (jl.ip . (jl.2p with the Allen-Cahn/Cahn- 
Hilliard equations. They have recently been studied in [9l fl~2| [IB]. 

Statement of the main result. Since we shall look for solutions, which are spherically 
symmetric, to problem (jl.ip - (ll.5p . the problem can be reduced to the one which is one 
space dimensional. To this end we now assume that the body force boundary and initial 
data and the unknowns, which are defined in the domain Q x (0, T e ), have the following 
form 

x 

b(t,x) = b(t,r)-, 7(£,x) = j(t,r), S (x) = S (r) 
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and 

u(t,x) = u(t,r) — , S(t, x) = S(t,r), 
r 

respectively, where T e is a positive constant which denotes the life-span of weak solutions, 
r = \x\, = {x € M 3 | a < r < d} for two positive constant a, d satisfying a < d, and 
6, 7, So are given functions and u, S are scalar functions to be determined, which depend 
only on t, r. We write 

x = (xj), u = (ui), T = (Tij), D = (Z>£j). 

Here and hereafter, i,j,k,l = 1,2,3, and we assume that D satisfies the properties of 
symmetry: 

0% = D% = D% = B%. (1.7) 
Moreover we assume that D satisfies 

Dfj = 0, if k + j; (1.8) 
Dl] = 0, if i + I, for any fixed j. Assume that D\- is independent of j; (1.9) 
Cu := Dfj. Assume that Cu is independent of I and is equal to \i. (1-10) 

3 

E kl := D pij = 0, if A: / /, and (1.11) 

Ekk is independent of k and is equal to A. (1-12) 



Under these assumptions equations (jl.ip - (|1 .3j) are reduced to 

d 2 A 29. 2 



+ = ^' (1-13) 

ar z r or r z 

d - f d 2 - \ d - 

^S+(-cu—S + T)\-S\ = 0. (1.14) 



dt \ dr 2 j dr 

Here J 7 , Q are nonlinear functions defined by 



g = g(^-s,b) = -^-s + -, (i.i5) 

Or) \ 

-A(— «+-«) +-De- eS + ${S) J , (1.16) 

[a t l^S,ls)=j^-—^-S. (i.i7) 

or or / r or 



The boundary and initial conditions become 



fi(t,r)=7(*,r), S(t,r) = 0, (t,r) E [0,T e ] x 90, (1.18) 
S(0,r) = 5 (r), reO, (1.19) 
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where ^(t, r) is defined by 7(i,r) = i(t,r)~. 

Remark 1. One can easily find an example which meets the above assumptions: The 
media is isotropic and homogenous. These assumptions still lead to an elliptic-parabolic 
coupled system, thus the reduced system possesses the main difficulties in the proof of 
existence of weak solutions as in Jlj/. 

To what follows, except Section 2 in which we reduce the problem to one dimensional 
form, we shall change the independent variable r to x, and drop the hat * for all quantities 
(except namely, u — >• u, b — > b, etc. Denote f x = gg/, fxx = J^/> e ^ c - The domain 
is reduced to an interval: Q = (a, d) is a bounded open interval with constants a < d. 
We write Qx e '■= (0, T e ) x $7, where T e is a positive constant, and define 

(v,ip)z = / v(y)ip(y) dy , 



z 

for Z = Q or Z = Qt £ - If u is a function defined on Qt b we denote the mapping 
x — > v(t,x) by v(t). If no confusion is possible we sometimes drop the argument t and 
write v = v(t). Since equation (|1.13p is linear, the inhomogeneous Dirichlet boundary 
condition for u can be reduced in the standard way to the homogeneous condition. For 
simplicity we thus assume that 

7 = 0. 

Then with these simplifications, equations (|1.13p - (|1.14p can be written in the form 

2 2 

u xx + —Ux = G, (1.20) 

x x A 

^S+(T-cuS xx )\S x \ = 0. (1.21) 
The boundary and initial conditions turn out to be 

it(t,z)=0, S{t,x) = 0, (t,x) G [0,T e ] x dQ, (1.22) 
5(0, x) = S (x), xen. (1.23) 

To define weak solutions of this initial-boundary value problem we note that because 
of \{\y\y)' = \y\ equation (jl.2ip is equivalent to 

^5-^(5,15,1^ + ^15,1=0. (1.24) 

Definition 1.1. Let b € L°°(0, T e , L 2 (VL)), S G L°°(Q). A function (it, 5) with 

M Gr(0,T e ;< M (n)), (1.25) 
5 e L°°(Q T J n L°°(0, T e , H^(n)), (1.26) 

is a weak solution to the problem \1.2U\) - (|1.23|) , if the equation M.2U\) is satisfied weakly 
and if for all (p G Cq°((—oo, T e ) x Q) 

(S,wh Te - <Px)q Te ~ (^|5,|, <p)q Tc + (S o ,<p(0))a = 0- (1.27) 



The main result of this article is 
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Theorem 1.1 To all So € Hq(Q) and b € C(Q T J with b t € C(Q T J there exists a weak 
solution (u, S) of problem hi. 20(1 - hl.23\) . which in addition to hi. 25]) - \1.27\ ) satisfies 



S t eLl(Q Te ), S x eLl(0,T e ;L°°(n)), (1.28) 

and 

(\S X \S X ) X €L*(Q T J, SxteLifaTjW- 1 '*^)). (1.29) 

Consequently we find spherically symmetric solution (u(t,r) -,S(t,r)) to the original 
problem (JT7E]) - ([T3]) . 

The remaining sections are devoted to the proof of Theorem 11.11 The difficulties 
in the proof of existence of weak solutions to the one dimensional problem are due to 
the following features: The system is of elliptic-parabolic type, it consists of a linear 
second order elliptic equation coupled with a nonlinear equation for the order parameter 
equation. The nonlinearity of this nonlinear equation is stronger than the one in [3J, 
where a one-dimensional initial boundary value problem is investigated. Moreover, this 
equation is degenerate and the nonlinearity depends non-smoothly on the gradient of 
unknown S. This can be judged easily from the fact that the coefficient v\S x \ of the 
highest order derivative S xx in the order parameter equation is not bounded away from 
zero and that it is not differentiable with respect to S x . 

The rest of this article is organized as follows: In Section 2, assuming that the domain 
CI, the elasticity tensor D and the misfit stain tensor satisfy suitable conditions, and 
that the solutions (u, S) to problem (jl.ip - ([1.5]) and the initial and boundary data are 
spherically symmetric, we reduce the original problem to the one dimensional form. 

Then to prove Theorem 1 1.1 1 we first consider in Section 3 a modified initial-boundary 
value problem which consists of (|1.20[) and the equation 

S t - cv\S x \ K S xx + T- (\S X \ K - k) = 0, x € Q, t>0 (1.30) 

with a constant k > 0. Here we use the notation 

\p\ K := ^Tp^. (1.31) 

Since (jl.30p is a uniformly parabolic equation we can use a standard theorem to con- 
clude that the modified initial-boundary value problem has a sufficiently smooth solution 
(u K ,S K ). For this solution we derive in Section 4 a-priori estimates that are uniform 
in k for k € (0, 1] . The assumption k £ (0, 1] is reasonable since we consider limits of 
approximate solutions for k — > 0. We shall see that the selection of a function in the 
form (11.311) results in a simpler proof of the existence of weak solutions than that in [3] . 

To select a subsequence converging to a solution for k — > we need a compactness 
result. However, our a-priori estimates are not strong enough to show that the sequence 
S x is compact; instead, we can only show that the sequence J Q Sx \y\dy = ^S£\S£\ has 
bounded derivatives, with respect to both x and t, in some suitable spaces, and thus 
can be proved to be compact. It turns out that this is enough to prove existence of 
a solution. For the compactness proof in Section 5 we use the Aubin-Lions Lemma; 
since one of our a-priori estimates for derivatives of the approximate solutions is only 
valid in L 1 (0, T e ; H~ 2 (ft)), we must use the generalized form of this lemma given by 
Roubfcek [24], which is valid in L 1 . 
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Despite we prove the existence of spherically symmetric solutions, the existence of 
weak solutions to the original problem (jl.ip - (jl.5p is still open. The method of the proof 
in this article and [3] (in which problem (II. ip - (II. 5p in one dimensional case is studied) 
is limited to one space dimension, since for the a-priori estimates it is crucial that the 
term | | in (ll.2ip can be written in the form ^(\S X \S X ) X . In the higher dimensional 
case the corresponding term \^7 X S\A X S cannot be rewritten in this way. Yet, we believe 
that these essentially one-dimensional existence results can also be helpful in an existence 
proof for higher space dimensions. 

2 Reduction to one dimensional problem 

We shall prove in this section that under suitable assumptions, the original problem can 
be reduced to a one dimensional problem. We now assume that the body force and the 
unknowns, which are defined in the domain Q x (0,T e ), have the following form 

b(t,x)=b(t,r)-, (2.1) 
r 

and 

u(t, x) = u(t, r)~, S(t, x) = S(t, r), (2.2) 
r 

respectively, where T e is a positive constant, r = \x\, £1 = {x € R 3 | a < r < d} for two 
positive constant a, d satisfying a < d, and u, S are scalar functions to be determined, 
which depend only on t, r, and b is a given function in t, r. 

Theorem 2.1 Suppose that the tensors D ande satisfy (jl.7p - (|1.12p . Then the following 
two statements are equivalent: 

1. (u, S)(t,x) of the form (|2.2p is a classical solution to the problem (jl.ip - (jl.5p with 
b chosen in (12. lj) . 

2. (u, S)(t,r) solves classically the problem (|1.13p - (|1.19D . 

Proof. To simplify notations, the Einstein summation convention applies to the rest of 
this section: When an index variable (e.g. k, I, but with an exception r in this article, 
for instance, S >r ^- in (|2.6p does not mean that we take the sum for the index r) appears 
twice in a single term that is a product of two or more numbers, it implies that we are 
summing over all of its possible values. However we shall still use the symbol S to avoid 
some possible confusion when in a single term, an index appears more than two times. 
Recalling 

x = (xi), u = (ui), T = (Tij), D = (D^). 

where i, j, k, I = 1, 2, 3. 

For partial derivatives, we denote for a function / = f(t, r) 

J,r — a i J,rr — o 9 • 

or or z 

An index j (or i, k,l,) after a comma in subscript of a quantity (for example, a function 
/ = f(t,x), a vector u = u(t,x) and tensor T = T(t,x), etc.) indicates the partial 
derivative with respect to Xj, namely 

df_ _ duj_ _ dT ik 

J.j — o j u i,j — o 5 ^-ik.i — o i ' ' ' • 
OXj OXj OXj 
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Similar convention applies to multiple indices after a comma in subscript of a quantity. 
We can thus write 



Xj, f Xi \ ^ijT Xj^Xj 



(2.3) 



rp . \ rp . rp . j\ . . rp^ rp . y . 

u—\ = n, r — 7T- + U— ^ -, (2.4) 

V / j T T 

XiXjXk „ T ((^Ej^Ej)^ $ij%k) 3XiXjXfc 

r ( 5ijXk (xjXy) fc) + 3xjXyXfc . . 

+u J - — , (2.5) 

Si = S r — . (2-6) 

' r 



Here 5ij is the Kronecker delta. 

Hence, the first two equations (jl.ip - (|1.2p can be rewritten as 

h Xl — n k K, _j_ n M -,, n fci F q 

_ * n fc* x i x j x k r / n fcl ( x i x j),k j^kl $ij x k _ n fcf SXjXjXk 

+^ (-r 2 (Z>g6ijx fc + DgCxiXj),*) + - D^S^. (2.7) 

From now on we take one £ from {1,2,3}. Invoking assumptions fjl .8|) - ()1.12j) we 
obtain 

3 

n fc< x i x j x k _ n jr7 XjXjXj Xi X\ X\ 

Since (xiXj)^ = S^Xj + Xi5jk, one has 

nfcZ ( x i x j),k jjkl 3ij x k _ jjkl SXjXjXk 



E/ SjjXj + ajj(5jj n jl $ij x j _ y.jl ^ x % x f 



3 



^ r 2 (3gj + 2gj) - 3x^r 2 _ 2r 2 ^ 



2 ^ 



8 



and 



r 



ij ^5 



3=1 
3 



E jl —r (25{jXj + Xi5jj) + 
ij ^5 



3=1 

E -r 2 (2xj + 3xi)+3r 2 Xj _ ~2xj 

i=l 

Using (|2.8p - ([2.10p . we are in a position to rewrite equation (|2.7p as 

b— = — ( fx(u trr H — & r ^ii) — XS r ) . (2-11) 



This holds, for r > a > 0, if and only if the following equation is satisfied 

/ 2 2 \ 
^t(n rr H — u r 9^)— XS r = b, (2-12) 

v ' r ' r J 

which is just (|1.13p . 

Next, to deal with the order parameter equation we make use of the following formula 
ip s (e,S) = -T-e + $'(S) = -De(Vu) ■ e + De ■ eS + $ {S). 
We evaluate De{Vu) ■ e. Invoking (|1.12p . 



De{Vu) ■ e = -D^u^e M + -D^u^eu 



j-, X{Xj -It? "ij tp XiXj 

u^ T hiij — k — h u Hi a hii ■ 



E * y u ,r-^-+u — - 3 — ) 



i=l 

= X (u r + , (2.13) 

where A is a constant as in (|1 . 12[) . Thus (|1 .3j) turns out to be 

S t + (-cvS )rr + T)\S )r \ = 0, (2.14) 

where T is the same function as in (|1.17p . Thus we obtain (|1.14p . 

To finish the reduction of the problem, we write the initial boundary conditions in the 
following form: j(t,x) = 7(t,r)^, Sq(x) = So(r). Thus we obtain the one dimensional 
problem and the proof of Theorem 12.11 is complete. 
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3 Existence of solutions to the modified problem 

In this section, we study the modified initial-boundary value problem and show that 
it has a Holder continuous classical solution, consequently we construct approximate 
solutions whose limit is a solution to the original problem (|1 .20j) - (jl.23p . To formulate 
this problem, let % € Cq°(M., [0, oo)) satisfy x(t)dt = 1. For k > 0, we set 



X«(t) := -X (- 

K \ K 



and for S G L°°(Q Te ,R) we define 

f Te 

(x K *S)(t,x)= Xn{t - s)S(s,x)ds. (3.1) 
Jo 

The modified initial-boundary value problem consists of the equations 

2 2 

u xx + ~u x = Q((xk* S) x ,b), (3.2) 

x x z 

St — cv\S x \ K S xx = —J--(\S x \ k — k), (3-3) 

which must hold in Qt s , and of the boundary and initial conditions 

u(t,x) = 0, S{t,x) = 0, (t,x) G (0,T e ) x dQ, (3.4) 
5(0, x) = S (x), xGO. (3.5) 

Now we want to rewrite the system as an equation with a nonlocal term. Applying 
the Sturm-Liouville theory for ordinary differential equations of the form 

-j- (p(x)y x (x)) + q(x)y(x) = 0, 
ax 

with suitable boundary conditions at x = a, d, we first solve u in terms of S x and b. To 
this end, we rewrite, by multiplying it by x 2 , (|2.12[) as 

L[u] := (jp{x)u x ) + q{x)u, (3.6) 
ax 

L[u]{x) = x 2 Q(( XK *S) x ,b), (3.7) 

and the boundary conditions are chosen as u(t, x) = at x = a, d. Here 

p(x) = x 2 , q(x) = —2. 

Consider the eigen-problem L[u] = au with a = and with u(t,x) = at x = a, d. 
It is easy to show from (|3.7p that 







[ L[u] ■ udx = I ( — (p(x)u x ) u + q{x)u 2 } dx = — / (x 2 u 2 + 2u 2 ) dx, 

Ja Ja \dx ) J a 



whence u = 0, and is not an eigenvalue of this operator. One asserts that for any fixed 
t G [0,T e ], there exists a unique solution u to ()3.2p . which can be represented by 

u(t, x) = f G(x, y) ( V -b{t, y) + -y 2 ( XK * S(t, y)) y ) dy. (3.8) 
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Here G(x, y) is the Green function, associated with the operator L, such that 

1. G(x, y) is continuous in x and y; 

2. For x^y, L[G(x,y)]=0; 

3. G(x, ■) = at x = a, d; 

4. Derivative jump: G'(y +0 ,y) - G'(y- ,y) = 

5. Symmetry: G(x,y) = G(y,x). 

Recalling the boundary condition (|3.4|) and integrating by parts we infer from 
that 

I rd \ f d 

u(t,x) = - G(x,y)y 2 b(t,y)dy / (G(x, y)y 2 ) y k * S(t, y)dy 

A* J a A* J a y 



i r d „, , 2 .,. , , a 



G{x, y)y z b(t, y)dy - - I 2G(x, y)yxn * S(t, y)dy 

A* J a A* J{x^y} 

— I G(x,y) y y 2 XK *S(t,y)dy. (3.9) 

Thus «(£, x) depends linearly on S and a nonlocal term of 5". 

To formulate an existence theorem for this problem we need some function spaces: For 
nonnegative integers m, n and a real number a £ (0, 1) we denote by C m+a (£l) the space 
of m— times differentiable functions on f2, whose m— th derivative is Holder continuous 
with exponent a. The space C a ' a ^ 2 (Q Te ) consists of all functions on Q T , which are 
Holder continuous in the parabolic distance 



d((t,x),(s,y)) := V\t~s\ + \x-y\ 2 . 

C m ' n (Q Te ) and C m+a,n+a / 2 (Q Te ) 1 respectively, are the spaces of functions, whose x- 
derivatives up to order m and i-derivatives up to order n belong to C(Q T J or to 
C a ' a ^ 2 (Q T J, respectively. 

Theorem 3.1 Let u, k > 0, T e > 0, suppose that the function b € C(Q Te ) has the 
derivative b t € C{Q Te ) and that the initial data Sq S C 2+a (Q) satisfy So\gn = So iX \gn = 
Soxx\dn = 0. Then there is a solution 



(u,S) G C 2 '\Q T J x C 2+a ^ 2 (Q. 



of the modified initial-boundary value problem h3. S\) - \3.5}) . This solution satisfies St x G 
L 2 {QT e ) and 

max \S\ < max |<So|- (3.10) 



Proof. Making use of (|3.9p . we rewrite the system (|3.2p - (|3.3p as a single equation 
St = ai(S x )S X x 

+a 2 [t,x,S,Sx,S, G(x,y)yS(t,y)dy, G(x,y) y y 2 S(t,y)dy (3.11) 

in Q Te , where S = Xk * S, 

ai(p) = cv\p\ K 
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and 



a 2 (t,x,S,p,r,s) = cT(t,x,S,p,r,si,s 2 )(\p\ K - re). 

Here J- is obtained by using formula (|3.9[) and inserting u, u x into the formula of T . 

Equation (|3.1ip is quasilinear, uniformly parabolic equation, which contains nonlocal 
terms. Then we can apply, with a little modification, [20^ Theorem 2.9, p. 23] to (|3.1ip to 
prove the existence of classical solution S K , and conclude that the estimate (|3.10p holds 
by applying the maximum principle to (|3.11|) . We refer the reader to the paper [3] for 
the details. Thus we complete the proof of Theorem 13.11 

4 A priori estimates 

This section is devoted to the derivation of a-priori estimates for solutions of the modified 
problem, which are uniform with respect to re £ (0, 1] . We remark that the estimates 
in Lemma 14.11 and Corollary 3.1, though stated in the one-dimensional case, can be 
generalized to higher space dimensions. 
In what follows we assume that 

0<re<l, (4.1) 

since we consider the limit re — y 0. The L 2 (0)-norm is denoted by || • ||, and the letter C 
stands for varies positive constants independent of re. We start by constructing a family 
of approximate solutions to the modified problem. To this end let T e be a fixed positive 
number and choose for every re a function Sq € C^°(0) such that 

\\S£ ~ S \\ H . {n) ^ 0, re^O, (4.2) 

where So € Hq(Q) are the initial data given in Theorem 1.1. We insert for So in (|3.5p the 
function Sq and choose for b in (j3.2[) the function given in Theorem 1.1. These functions 
satisfy the assumptions of Theorem 13. 1\ hence there exists a solution (u K ,S K ) of the 
modified problem (|3.2p - (|3.5p . which exists in Qt e - The inequality (|3. 1Q[) and Sobolev's 
embedding theorem yield for this solution 

sup \\S K \\ L oc( Q ) < sup \\Sg ||l<x.(q) < C. (4.3) 

0<K<1 0<K<1 

Remembering the formula (|3.9p and assumptions of b, we show easily that u K be- 
longs to C 1,1 (QT e )) and conclude from (I3.2p that also ||u£||£<x>(q t \ < C, and invoke the 
definition of T\ to get 

max|J-i(u K ,<,S K )| < C. (4.4) 

With the help of this estimate, we can evaluate derivatives of S K . 
Lemma 4.1 There holds for any t € [0, T e ] 

||S^)|| 2 + c^ f [ \S K x \ K \S* x \ 2 dxdr<C. (4.5) 
J o Jn 
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Proof. From the assertion that S^ x € L 2 (Qj> e ) in Theorem 13. 11 it follows that there holds 
for almost all t 

~\\S£(t)\\* = j^(t)S K xt (t)dx. 

Making use of this relation and (|4.4p we obtain by multiplication of (|3.3p by —S xx and 
integration by parts with respect to x, where we take the boundary condition (|3.4p into 
account, that for almost all t 

^^ll^a-ll 2 + J ci, \^x\k.\^xx\ 2( ^ x = J J~(\Sx\k ~ K )S xx dx 

< C f (l + \S«\)(\SZ\ K + K)\S% x \dx 
Jn 

= C[ I \S X \ K \S xx \dx + / K\S xx \dx + / \S X \ \S x \ K \S xx \dx / k\S x \ \S xx \dx J 
\Jn Jn Jn Jn J 

= C{h+h + h + h) (4.6) 
Now we estimate Ii (i = 1,2,3,4). For I±, we have 

f K | K \ 

Jn 

< f I \SZ\ K \SZ x \ 2 dx + C u [ (\S:\ K fdx+C, (4.7) 
° Jn Jn 

where we denote C v a constant depending on v. By definition, there holds 

\S K X \ K > k. (4.8) 

Thus we can use the second term on the left hand side of (|4.6p to absorb I<i- By the 
Cauchy-Schwarz and Young inequalities, 

h < C [ K^(J\S^ x \)dx <c( [ Kdx] 2 ( [ K\S^ x \ 2 dx) ' 
Jn \Jn J \Jn J 

< ^ ! \S« x \ 2 dx + C v . (4.9) 



XX I 

n 



Here we have used the fact that < a < x < d which implies the term ^ contained in T 
is uniformly bounded from below and above, and C U K is a constant depending on u, k. 
Moreover, /3 is evaluated by 



i i 



h < C j^\SZ\\SZ x \dx<C^K\S«\ 2 d x y (J K \S^\ 2 dxY. (4.10) 
Using the Nirenberg inequality in the following form 

H/zll < C||/ asB ||i||/||f 00(n) + C7'||/|U» (n)j (4.11) 
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one infers from (I4.10P that 

< c K *(\\s2 x \\s + i) ( [ K\s* x \ 2 dx 



i 

3 i n'\\ Q K \\ 1 / / *-l Q K \ 2 r 



= Ck(\\S^ +1 + \\S^ 
< ^J^ x \ 2 dx + C u , (4.12) 

4 . 

where we used the Young inequality of the form ab < 5a 3 + C$b for non-negative real 
numbers a, b. Here Cg is a constant depending on 5. 

I3 is the most difficult term to deal with. Again by the Cauchy-Schwarz inequality 
one gets 

h = C I W\\S^\l){\S^\\\S^\)dx 
Jn 

< c Qf \s%\ 2 \sz\ K dx^ 2 |^u|5-j 2 d^ 2 

= CJ§ • (7 |S£U^J 2 ^V (4.13) 

To deal with J, we recall the boundary conditions for S K and rewrite by integration by 
parts 

J = I {S" x ) 2 \S*\ R dx = [ S«(S«\S«\ K )dx 
Jn Jn 

S K S« x {\SZ\ K + S«(\y\ K y\ y= s s )dx 
= — f S K S XX ( \S X \ K + — J 

(4.14) 

Applying estimate (|4.3|) and invoking the definition of from (|4.14p one obtains 

J < C [ 2\S K x \ K \S K xx \dx 
Jn 



< C( I \S K x \ K \S^ x \ 2 dx) [ I \S%\ K dx 
n 



1 1 

2 / f \ 2 



< CU\S^\ K \S^dx) + (4.15) 
Therefore, (|4.13p becomes 

h < C(||^|| + l)l ( I \S K x \ K \S K xx \ 2 dx 

\Jn / 

< C ^ JjS*\ K \S*J 2 dx + C u (\\S K x \\ 2 + l). (4.16) 



-+- 
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Here we have again used the Young inequality: ab < 5az + C^b A . 

Combining estimates (|4.7p - (|4.16p . subtracting the terms ^ f n \S x \ K \S xx \ 2 dx and 
f n \S xx \ 2 dx on both sides of inequality (|4.6p . splitting the second term on the left 

hand side of (|4.6p into two equal terms, and recalling the property (|4.8p . we derive 



IA||QH|2 + ^ 
2dt" x " 4 



I OKI I C/t |2 J , 

Pi « Past ux ' 



CZ/K 



>S^,^, I dx <C C 1 1 5'^, 1 1 ^ ~t~ Cy . 



(4.17) 



Then using the Gronwall inequality, one gets (|4.5p . noting also (|4.2p . And the proof of 
this lemma is thus complete. 

Furthermore, we obtain 
Corollary 4.1 There holds for any t € [0, T e ] 



o Jn 

t 



o Jn 



(\SZ\ K \SZ X \)3 dxdr < C, 
{\S K x S K xx \)Uxdr < C, 



\y\ndy 







dr < C, 



\y\ndy 



dr < C, 



L°°(f2) 



I OK I OK II 

I XI X H L3 ( 0]Te . L oo(Q)) 



I C«|| 3 

i J i iiL°°(n) 



< c, 

dr < C. 



Proof. For some 2 > p > 1 we choose q, q' such that 

2 1 1 
?=- - + - = 1. 

By Holder's inequality, we have 

f / (\S*\k\S K xx \T dxdT 

Jo Jn 

(i^u) f (d^u) § i^r) 



dxdr 



< 



< 



o Jn 

t 



o JJ) 



o Jn 



(|S£U)^ dxdrY U J \S^rdxdr\ 
(I^U)^d^r) 2 ^ j^\ K \S^ x \ 2 dxdry . 



(4.18) 
(4.19) 

(4.20) 

(4.21) 
(4.22) 

(4.23) 



(4.24) 



Estimate ()4.5p implies that if p satisfies < 2 (i.e. p < |) then the right hand side of 
(03]) is bounded. This yields estimate (ET181) . Then (1439]) follows from ()4TT8D and the 
estimate (|4.5|) . and the fact that \S X \ < \S X \ K . 
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Next we are going to prove (|4,20p . Writing 



I OK I QK 
\&X \K D XX 



\y\ndy , 







(4.25) 



and invoking that the primitive of \y\ K is equal to 



- iy\Jy 2 + k 2 + k 2 log (y + V y 2 + k 2 



which, thanks to logx < x — 1 for all x > 0, is bounded by C(y 2 + 1), we then show 
easily that 

\ K dydx<C / {\S*\ 2 + l)dx < C. 
Jn 



n Jo 



To apply the Poincare inequality of the form 

11/ - f\\LP(Q) < C\\fx\\lA>(Ci) 



where / := m J n f(x)dx, we choose 



P=3, / 



\y\ K dy, 



and obtain 



\y\ndy 



4 
3 

L3(Q) 



dr 



c 



\y\ K dy 



4 
3 

L3(fi) 



dr + C 



r/r 



L3(f2) 



< C / \\\S£\ K S^\\\ dr + C 



(4.26) 



which implies, by (|4.18p . that 



\ K dy 



4 

3 

L3(Q) 



dr < C. 



(4.27) 



Hence (EL2U1) follows, and we get f Q S * \y\ K dy € Li (0, T e ; VF 1 '! (0)). Making use of the 
Sobolev embedding theorem, we get (|4.2ip . 

It remains to prove estimate (|4.22p . since (|4.23p is equivalent to ()4.22p . We rewrite 
Jo x \y\ndy as 



si 



\ndy 



\y\dy + 



s';. 



\y\y 





+ 





U - \y\)dy 

K- 2 

dy 



o li/U + M 



\S£\sz + 



o \y\n + \y\ 



dy. 



(4.28) 
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Thus 



\vUdy 

/ x 



Ore i I Ck\ 



(4.29) 



By 
that 



and the Young inequality we obtain from (j4.5j) and the assumption that k < 1 



k2 $xx 



\Ck\ _i_ I OkI 
\ v k ' .t 



< Ik^^jI, thus 



< I J (k 2 + k\S^ x \ 2 ) dxdrj < C. 



(4.30) 



Combination with (f4T20|) and (f4729|) yields 



\S^ \S^) X \\ 4 , „ < c 

1 xl x; U LS(Q Te ) - 



< C. (4.31) 



It is clear that \S%\S% < C J n \S^\ 2 dx < C. Applying again the Poincare inequality 
to the function / = \S X \S X , we arrive at 



Hence this, combined with (|4.3ip . implies that 



| S 1 ^ | tS"^ | 



4 ,4 < C, 



one concludes by using the Sobolev embedding theorem that 



which is 



I OK I OK II <C C 1 

1 xl x "L*(0,T K ;L°°(O)) ~ ' 



1 x "L3(0,T e ;L°°(n)) ~ 



This completes the proof of the corollary. 

To apply some compactness lemma to the approximate solutions, we need estimates 
on the time derivative of the unknown S K and also \S£\S£. 

4 

Lemma 4.2 The function belongs to Ls[Q Te ) and we have the estimates 



\\ S t\\L*/Z(Q Te ) ^ C ) 
Pxl'5 r a)tlll,i(o,T e ;fl'- a (n)) - ^ 1 



(4.32) 
(4.33) 



Proof. From equation (|3.3|) and the estimates (|4.18p . and (|4.5p we immediately see 

4 

that € Lz(QT e ) and that (|4.32p holds. Therefore we only need to prove the second 
estimate. 
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To verify (|4.33l) we must show that there exists a constant C, which is independent 
of k, such that 



((\ S x\ S x)t^)Q T ^ C \\ i P\\L°°(a,T e ;H2(n)) 



(4.34) 



for all functions 99 £ L°°(0, T e ; Hq To prove ()4.34p . we first prove that for any 

1 > 5 > there holds 



< C'||(/?|| L oo( 0]Te . H 2(^)) 



(4.35) 



for all functions € I/°°(0, T e ; //q (il)). Here <5 is independent of k. Inequality (|4.34[) is 
obtained from this estimate as follows: From S£ € L°°(0, T e , L 2 (0)) C L 2 (Q Te ), S xt € 
L 2 {Qt e ) and | \y\$ — \y\ \ < 5 — > as 5 — > we infer that || \S£\s — \S£\ \\io C ^q t ) — >• 0. A 
straightforward computation yields that 



Is/Ms/ 



I QK I QK 

l°X \S D xt 



(4.36) 



Therefore, ( Jjf* \y\&dy) = \S£\ s S£ t -)• S£ t strongly in L 2 (Q Te ). Whence, as 5 0, 



for all ^ e L°°(0,r e ;i^(n)) C L°°(Q T J. This relation together with ([435]) implies 
(Oil) . 

Thus it suffices to prove (14.35p . To simplify the notations we define 



TZ K :— cv\S x \ K S xx J- K (\S^\ K k). 
Here T K = T{u R , tt£, S K ,S£). Recalling estimate (|4.4p . we have 

< C(\S^\ K \S^\ + (1 + |S£|)(|S£U + nj) . 



(4.37) 



(4.38) 



Multiplying equation (|3.3p by (\S£\sip) x , integrating the resulting equation with respect 
to (t, x) over Qx e , using integration by parts for the term with the time derivative and 
noting (|4.36p . we obtain 







(S t K -K K ,(\SZ\ s ip) x ) QTe 

~ i S xt, \ S x\d<f)Q Te ~ {R-k, {\Sx\s<p) x )q t 



\y\sdy) ,<p 



n K , (\y\5)'\ y=s . S^tp) - (K K , \S2\ s ip x ) .(4.39) 



Remembering that S xt € L 2 (Qt e ) for any fixed k, we see that the first term in the second 
equality of (|4.39p is properly defined. 

To estimate the last two terms on the right hand side of inequality (|4.39p . we note 
that there holds 

y 



y\s 



\y\s 



< 1 and \y\ s < \y\ + 1, 
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which yields the estimates 

K K , (\y\s)' 



< (\n K \,\s^\) 



and 



< (|s*U|s*J 2 , M) QTe + ((1 + \SM\SZ\k + «), \s^<p\) QTe 

< C||(^|| L o (Q T J + I < C\\(f\\ Lac{0jTe . H 2 {n )) + I, 

{n K ,\S2\ s <p x ) QTe < c[ (\S^\ + l)\S^\ K \S^ x \d(r,x) 

+ C [ (l + |^|) 2 (|S£| K + K ) \ip x \d(r,x) 



>Q Te 
Jl + J 2 - 



We estimate I first. Write 



I = c(\sz\ k + k+\sz\\sz\ k + k \sz\,\sz x <p\) Qt<: 
= h + h + h + h. 

One needs to estimate Ii (i = 1,2,3,4). By estimate ()4.5h . I\ can be treated as 
h = C(\S«\ K ,\SZM) QTe 

< C [ \SmS K x \l\S K xx \W\d{t,x) 

JQ Te 

< C r || |^|||| £ 4 (n) || \S^J \\ V \\ LHQ) dr 

Jo 

i i 

< c(J "lll^ll^H 2 ^" 

< C||^llL 2 (0,T e ;L 4 (n))) 

and 



llvlli*^)^ 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



h = c(K,\s«M) QTe 

< C [ K \SZ x \\<p\d(t,x) 

< CK\\S^ x \\ L 2 { Q Te) \\ip\\ L 2 { Q Te) 

< C\\<p\\ L 2 {0tTe . L 2 m . (4.44) 
With the help of (|4.19p and of the Cauchy-Schwarz inequality, we deal with I4 as follows 

h = C (k\S«\,\S«M)q t 



< C £ U\SZS«Jdx 



\^dx dt 



< C\\S^S^ x \\ L 4 {Q ^ )m \ L 4 {QTe) 

< C|MU 4 (0,T e ;L4(n))- 



(4.45) 
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The remaining term I3 is the most difficult to evaluate. Making use of estimates ([4.5 
and (|4.23p . we have 



h = c(\s£\\S2\ K ,\s2 x <p\) 



< C\\<P\\L°°(Q Te ) 



< C||^|| L oo ( Q 2 



\Jn 



( / \sz\ K \sz x \ 2 dx) dt 



(ll^|II,» (n) + i)( / \S K x \ 2 dx 



'(iraicc (n) + i) 2 <a 



\s^\ K \s^\ 2 dx) dt 



Jn 



< C\\<Ph°°{Q Te ) 

< C|MIl-(q t j- 

Next, we consider J±, J2. The term J\ can be handled as 

Ji = C [ (\S*\ + l)\S«\ K \S^ x \d(T,x) 



\S^\ K \S^ x \ 2 dxdtj dt 



(4.46) 



< c||^|| £0 c WTe) j o U(i + \sz\ydx 

< c\\<p x \\ L oo(Q Te ) U + \\s x \\ 2 Lac{n) \\s x \ 







(\SZ\ K \SZ x \)sdx) dr 
{\S K x \ K \S« x \)*dx\ dr 
(\S«\ K \SZ x \)idxdr) 







< C||v3 X .|| L 00(Q Te) 

< ClI^IU-ccTe;^^)). (4.47) 

Here we used the estimates in (j4.5|) and Corollary 14.11 which will also be used to evaluate 
the term J2. Invoking inequality (j4.38[) . we obtain that 

J2 < C [ (1 + |S£|) 2 (|S£U + k) \<p x \d(r,x) 



< C||^|| Loo(QTe) / (l + \S^)d(r,x) 

JQ Te 



< C II <Px II L°°(Q Tr ) 



1 + ||5, 



\S^\ 2 dx ) dr 



< C\\ ( Px\\L°°(Q Te ) \1 + /_ ll^||i°°(Q)^ 

< C||99|| iO o( 0iTe;Jf 2( n )) 



(4.48) 



Combination of (|4.39p - (|4.48[) and using the Sovolev embedding theorem yield 

\y\sdy) ,<p 







< C {\W\\l°°{0,Tz,HI{Q)) + IM I L°°(Q Te ) + ll^llL4(0,T e ;L4(C))) 

< cikll 



(4.49) 



which implies (J4J35J) and we complete the proof. 
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5 Existence of solutions to the phase field model 

We shall use in this section the a priori estimates established in the previous section to 
study the convergence of (u K , S K ) as k — > 0. We shall show that there is a subsequence, 
which converges to a weak solution of the initial-boundary value problem (|1.20p - () 1 . 2 3 1) , 
thereby proving Theorem ll.il 

Note first that the estimates in Corollary 14, 1 1 and Lemma [4. 2\ the fact that f2 is bounded, 
and Poincare's inequality imply 

ll^ K |lw 1 . 4 / 3 (Q Te ) - C i (5.1) 

for a constant C independent of k. Hence, we can select a sequence n n — > and a function 
S e W^iQrJ, such that the sequence S Kn , which we again denote by S K , satisfies 

~~ ^IIl 4 / 3 (Qt £ ) °' S£ — v S x , S£ — ^ S t , (5.2) 

where the weak convergence is in L 4 / 3 ((5t e ) . 

As usual, since equation (|3.3p is nonlinear, the weak convergence of S£ is not enough 
to prove that the limit function solves this equation. In the following lemma we therefore 
show that S£ converges pointwise almost everywhere: 

Lemma 5.1 There exists a subsequence of S£, we still denote it by S£, such that 

a.e. in Q Te , (5.3) 
a.e. in Q Te , (5.4) 

4 

weakly in La(Q T J, (5-5) 
strongly in lI (0, T e ; L 2 (ft)), (5.6) 

strongly in L% (0, T e ; L 2 (n)), (5.7) 
as k — y 0. 

The proof is based on the following two results: 

Theorem 5.1 Let Bq be a normed linear space imbedded compactly into another normed 
linear space B which is continuously imbedded into a Hausdorff locally convex space 
B\. Assume that 1 < p < +oo, that v,Vi € L p (0,T e ; Bq) for all i € IN, that the se- 
quence {vijieiN converges weakly to v in L p (0,T e ; Bq) and that {^}j 6 ]N is bounded in 
L l (0, T e ; B\). Then vi converges to v strongly in L p (0,T e ; B). 

Lemma 5.2 Let (0,T e ) x £1 be an open set in M + x M. n and assume that 1 < q < oo. 
Suppose that the functions g n ,g € L 9 ((0, T e ) x Q) satisfy 

\\9n\\Li((o,T e )xn) <C, g n ^r g almost everywhere in (0,T e ) x O. 
Then g n converges to g weakly in L q ((0,T e ) x Q). 



I <? K I — \ I <? I 

I Q K I _i. I Q I 
\ u x \ K I x I ' 

\y\dy ->■ 2^ 1*9x1, 



\y\ndy 



2^x\S x 
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Theorem 15.11 is a general version of Aubin-Lions lemma valid under the weak assumption 
dtVi € L 1 (0, T e ; B\). This version, which we need here, is proved in Simon [25] and in 
Roubfcek [24] . A proof of Lemma 5.2 can be found in [214 P-12]. 

Proof of Lemma 15. j I We choose p = | and 

B = W 1 ^(Cl), B = L 2 (Cl), B l = H~ 2 (Cl). 

These spaces satisfy the assumptions of the theorem. Since the estimates (I4.18P , ([4,20j) 
and (|4.33p imply that the sequence J Q X \y\dy is uniformly bounded in L p (0,T e ; By) for 

k — y and (^Jq x \y\dy^ is uniformly bounded in L 1 (0, T e ; B\), it follows from Theo- 
rem 3.1 that there is a subsequence, still denoted by J Q x \y\dy, which converges strongly 
in L p (0,T e ;B) = Ll(0,T e ;L 2 (Cl)) to a limit function G € lI (0, T e ; L 2 (Cl)). Conse- 
quently, from this sequence we can select another subsequence, denoted in the same 
way, which converges almost everywhere in Qt £ - Using that the mapping y \— > f{y) := 
Jo \£\d£ = \y\v\ nas a continuous inverse f^ 1 : R — )■ R, we infer that also the sequence 
= / _1 ^ |y|<iy^ converges pointwisely almost everywhere to f~ 1 (G) in Qx e - From 
the uniqueness of the weak limit we conclude that f~ 1 (G) = S x almost everywhere in 
Qt £ - 

For the proof of (|5.7p we write 

r s x r s ': r s x 

/ \y\ndy = \y\dy + (\y\ K - \y\)dy = h + h- 
Jo Jo Jo 

It is easy to estimate I 2 as ||/ 2 ||l 2 (q t J < Wl 2 (q Te ) < Ck\\S% \\ L ™(o,T e ;L 2 (n)) <Ck^ 
0. Therefore, J Q X \y\ K dy — > lim^oA = 5 1 <S"x | •S'x strongly in La (0, T e ; L 2 (Cl). This is 

To prove (|5.5p we note that the estimate \S£\ K < \S£\ + k and the inequality ()5.ip 

4 

together imply that the sequence \S^\ K is uniformly bounded in La(Q Te ). Thus, (|5.5p is 
a consequence of (|5.4p and Lemma 5.2. 

Proof of Theorem 1.1: Define the function u by inserting S into (|3.9|) where S 1 is the limit 
function of the sequence S K . We shall prove that (u,S) is a weak solution of problem 

(oun - dm}. 

Remember first that by Lemma 14.11 we have S € L°°(Qy e ). From this relation, from 
the above definition of u we immediately see that u satisfies (jl.25|) . Observe next that 
\\^ K \\L ao (o,T e ;H^(Q)) — by Lemma I4TT1 and Sobolev's embedding theorem. This implies 
S G L°°(0, T e ; Hq (CI)), since we can select a subsequence of S" 1 which converges weakly 
to S in this space. Thus, S satisfies (jl.26p . 

Noting that from ([3~Tjl and (Oil 



< ll(s-x«»s)y (QrJ + ||s«-siy WTi) ^o, (5.8) 

for k — )■ 0, we conclude easily that the function u defined in this way satisfy wealkly 
equation (|1.20p . It is thus enough to prove that the equation (|1.20p - (jl.2ip are fulfilled 
in the weak sense. By definition, these equation are satisfied in the weak sense if the 
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relation (ll.27p holds. To verify (II. 27ft we use that by construction S K solves ()3.3p . Now 
we multiply equation (|3.3p by a test function 93 € C^°((— oo,T e ) x f2) and integrate the 
resulting equation over Qr e , then obtain 

= (S?,ip) QTa + (-cv\SZ\ K SZ x + F K (\SZ\ K -K),ip) QT 
= -( s o,<f( Q )h- (S K ,(ft)Q Te + lev \y\ K dy,(p x ) 

\ JO J Q Te 

+ {F*{\S*\ k -k)^) Qt . 
Equation (I1.27P follows from this relation if we show that 

(S K ,V?(0))n -> (5o,<^(0))n, (5.9) 



{S K ^ t ) QTe (S,<p t ) QTe , (5.10) 



y\Kdy,(fx) -> (~\S x \S x ,<p x \ , (5.11) 
{F k {\S«\ k -k\v) Qt -> (T,<p) Qt , (5.12) 







for k — y 0. Now, the relation (|5.9[) follows from (|4.2p . and the relation (|5.10p is a 
consequence of (|5.2p . By (|5.7p . one has ()5.1ip . 

To verify (pU2j) we note that (pT8]) . P~23|) and the definition of J" K yield 

ll^(I^L-*)ll Lf(QTj < C, (5.13) 

F K (\S£\k — «) -> ^ I <$e I, almost everywhere. (5-14) 

Then by Lemma 15.21 

4 

weakly in L3(Qt £ ), which implies (15.12|) . Consequently (jl.27p holds. 

It remains to prove that the solution has the regularity properties stated in (|1.28p 
and (|1.29p . The relation St € Ls(Qx e ) is implied by (|5.2p . To verify the second assertion 
in (jl.28p . we use estimate (I4.23P to get 

<"T e 8 

II <? K II 3 rlt < C 

\\°x\\L°°(n) u,L — u - 

This inequality and S x in L§ (0, T e ; imply S x € L§ (0, T e ; 

To prove (|1.29p . we recall that Jq 5 * |y| K <iy converges to |<S , X |S' X strongly in the space 

Ls (0, T e ; L 2 (S1)) C L3(Q T J and that ( J Q 1 |y| K <iy) is uniformly bounded in Lz (Qr e ) 

x 4 

for k — y 0, by (|4.18p . This together implies that (\S X \S X ) X € L3(Q Te ). Finally, to prove 
the second assertion of (jl.29p we choose a test function ip € L 4 (0,T e ,W ' multiply 
equation (13. 3p by — and integrate the resulting equation over Qt s to obtain 

= (S? - K K , -<p x ) QTe = (S K xt , v) QTe + {n K , <p x ) QTe , (5.15) 

with 1Z K defined in (|4.37|) . Invoking the estimates (|4.5p and (|4.23p we deduce that 

11 ""(Q Te ) - ' 
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hence equation (|5.15p yields 



{Sxt^) QTe < ll^ll L | (QTe) II^IU 4 (QT e ) ^ c y\\mo,T e] w^m • 

4 -.4 

and this means that S£ t is uniformly bounded in La (0, T e ; 14 7- '3 (f2)). From this estimate 

4 

and from — St in L^(QT e ) we deduce easily that 5 x t belongs to the dual space of 
L A (0,T e ; Wq A {Q)), which is Ls (0, T e ; VF" 1 ^ (fi)). 
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